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Oblique propagation of electromagnetic waves in an inhomogeneous magneto-ionic medium is 
described by a pair of coupled, linear, second-order differential equations, with coefficient functions 
having one essential singularity and more than one pole. An analytic method to study this general 
type of propagation has been developed in this paper. 

Mathematical expressions for all the physical observables in the ionospheric investigations have 
been derived for a parabolic electron layer. The expressions for the phase retardation and the tilt 
of the polarization ellipse differ radically from those derived from the ray optics of Appleton and 
Hartree. 

1 . Wave Equations 

The general type of magneto-ionic propagation is described by the equation 



£+* + «„ x 



i) 2 E 



m dt 



(i) 



where mks units and symbols (see appendix) have been used. 

The statement of the physical problem and a brief summary of the previous work in this line 
have been discussed in part I of this paper. 

In the above wave equation let, (i) the partial time derivatives be replaced by ia>; (ii) the 
steady magnetic field H ( >, (coh = e/JLoH e /m) be in the xz-plane with the vertical z-axis as the direc- 
tion of propagation; (iii) the collision frequency v be independent of z; and (iv) /V, the electron 
density, be given by 



N(z) = Njl~) for z*<P, (2) 

= for z 2 > 



where z = is fixed at the layer of maximum ionization. 

Under the above set of assumptions (1) decomposes into 



^+[Ki(z)E x -iQUz)Ey]=0, ^+[K&)E y +iQUz)E x ]=0, (3) 
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Z=\ — ivla), $=— + 8, 8 = magnetic dip angle, 
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For = — (transverse propagation), the above pair of coupled equations reduces to a pair of 

separate equations, discussed in part I of this paper. In the longitudinal case with = or 77% 
the linear combinations (E x ±iE y ) obey separate equations similar to the one satisfied by the trans- 
verse ordinary or o-mode. Except for these limiting cases, separate wave equations obeyed by 
linear combinations of the type (f\{z)E x +f 2 (z)E y ) cannot be set up for any type of functional de- 
pendence of N with z. Hence general magneto-ionic propagation is always governed by coupled 
wave equations. 

2. Nature of the Solutions 

For frequencies used in radio communications co# <^ co, i.e., \z 2 \ <^ l 2 y and hence solutions for 
the region z 2 > \z&\ are of more practical importance than those for the magneto-ionic core with 
z 2 < \zl\. Changes in the variables from E x , iE y , and z to ai = £ 1/4 ; E x , u 2 = £; 1,4 E y , and g=iKz 2 
reduce (3) to the form 
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If U\—* uioo, 112^ U2cc as £~ >0 °, the equations obeyed by the asymptotic functions Ui^ U2 X , 
when the terms of the order of l/£ 2 and smaller are neglected, assume the form 



#1*100 / 1 2fi\(Pu laa n-Q 2 n 
<%* \ 2 t) de L 16 2f 



^100 = 0, 
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Solutions of (5) can be expressed as 



Wlo 



= / c e ^n^~^ )aj ^' 



(6) 
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where 



aj = -(l+n/2r,j), r h , 2 = ±hVT+Q, 



Vs 



4=± i jVr^, 



and the contour C is a combination of four loops, C/s,j= 1 to 4, each round one of the four singu- 
larities Tjj = ±Vl ±Q of the integrand. Along any ray in the 17-plane, the contribution to (6) from 
the loop Cj can be written as 



Jc J 



nV 



■ VK) a K 
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where P\j is the phase factor independent off. Evaluation of these phase factors from the asymp- 
totic solutions of the type (6) is necessary for the study of the Stoke's phenomena and the esti- 
mation of the observables. The solutions of (4) can now be expressed as 



ut = 2 Puem* Fij(Z ), "2 = X iVV F2M ). 



(6a) 



where Fij(g) and F 2 j(f ) are series in descending powers of f with unity as the leading term. These 
series can be evaluated by substituting (6a) in (4) and equating the coefficients of the different 
powers of (1/f) for each of the four wave terms e 17 /, after l/(f — f ) has been expanded. Compari- 
son of these solutions with th ose fo r the transverse modes reveals that the wave factors exp 
[+ i Vl + ^f] and exp [— Jvl + ^f] represent respectively the upgoing and downcoming com- 
ponents of the o-mode while exp [Jvl — Q £] and exp [— 2 Vl _ Q £] represent those components 
of the e-mode. 

The four observables, viz, R, the ratio of amplitude of the downcoming to that of the upgoing 
wave, 0, their relative phase difference, ( M, the ratio of the axes of the polarization ellipse of the 
downcoming wave, and ^, the tilt of the major axis of that ellipse are related to these solutions as 



fte-** = 



PiiFutf) 
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z =-i for the o-mode, 



for the e-mode, 



&fc- { *-- 
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PllFu z=-l 

and V= i tan" 1 [2^ sin iff 1(1 -i^ 2 )] 



for the downcoming o-mode, 



for the downcoming e-mode, 



(7a) 
(7b) 
(8a) 

(8b) 
(8c) 



Evaluations of these solutions in (6a) and expressions in (7, 8) complete the study of the general 
magneto-ionic propagation in a parabolic ion layer. Before passing on to the evaluations of the 
expressions for the observables, the following salient features of the propagation as predicted by 

the above analysis need to be stressed: 

(i) There ar e only two modes of propagation characterized by the indices tj^ 2 =± 2 Vl + () 
and T7 3 4 =± \ Vl — (), with the two signs of each mode giving the waves traveling in opposite 
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directions for that mode. This result holds for any type of functional dependence of TV and v with 
z. Hence the presence of a third echo (triple-splitting observed in polar latitudes) is ruled out 
from the general magneto-ionic propagation. 

(ii) The complex polarization ratio &te~^ for the upgoing wave of either mode is not the same 
as that of the downcoming wave of the same mode. If it were so, suitable linear combinations of 
the dependent variables would give two decoupled second-order wave equations. 

(iii) The sharing of the energy between the two modes is controlled by the parameters of the 
problem and the propagation angle. This gives rise to a new observable. 

3. Expressions for the Observables 

The nature of the analytic solutions being known, a method based on the basic characteristics 
of (3) and their solutions (6a) will now be developed for the evaluations of the expressions for the 
observables. 

Writing iE y — p{z)E x in (3) and eliminating E x and iE y from the pair of equations by raising 
the order of differentiation, it can be shown that p(z) satisfies the equation, 

2pp'"-3p M2 + 4p' , [p(K 1 -^ 2 ) + (?L(l-p 2 )] 

- 2p' 2 [3K 1 -K 2 - 4p<?L] - 2p'[p(K{ - KQ + QL'(\ - p 2 )] 

- [(£, - K 2 )p + QL(1- p 2 )] 2 = 0, (9) 

where the dashes represent the derivatives with respect to z. 

Corresponding to each solution pj of (9), the components E x and E y of the electric vector sat- 
isfy the first-order equations 



dE x 
dz 



>;-{(*, -* 2 ) Pj +(?L(i-pj)} - 

2pj 



E x = 0, (10) 



and 



dE y 



p'j-iO^-Kjpj + QLa-pj)} 2pj 
2pj Pj 



dz 

Replacing iE y by pE x , the first of the equations in (3) can be written as 

d*E x 



dz 2 



?+[Ki-pQL]E x = 0. (11) 



Since for parabolic layers, solutions of (3) are symmetrical about z = and have an essential 
singularity of grade 2 at °° together with a pair of regular singularities at z 2 = zl, the expression 
(Ki — QpL) in the above equation must fulfill the following conditions: 
(i) It must be symmetric about z = 0; 

(ii) It must not contain powers of z higher than the second; 
(iii) It must not have any pole other than those at z 2 = z%. 
From the expressions of K\ and L, it follows immediately that p(z) must be of the form 

Substituting (12) in (9) and equating the coefficients of different powers of z 2 and (z 2 — z 2 ,) 
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to zero, it can be shown that 



1, fa = h 2 h 3 l 2 sin 2 0/2 cos = 



/WO)// 



V 



sin 2 6 



\ w 2 „ / \2 cos d 



(13) 
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1 



sin 2 



and for all values of n > 2, /3 n 's are single-valued. 

Since for each root of q, there are two roots of fa and for each pair of values of q and yS 2 all 
the other /3 w 's are single-valued, (9) has four and only four distinct solutions. Corresponding to 
each solution of (9), a pair of solutions for E x and E u can be evaluated from (10). The two roots 
of q characterize the two distinct modes of propagation, while for any of these modes, the two roots 
of fa specify the upgoing and downcoming waves bejonging to this mode. The method outlined 
here is exact and not approximate, though it is not applicable for values of the parameters h 3 and 
6 which make the series (12) divergent. For ionospheric propagation the method is applicable 
to all latitudes excepting a small belt of 5° around the equator. 

The complex phase factor Re~ i(i) between the upgoing and downcoming waves at the bottom 
of the layer depends on the asymptotic solutions of (11). For such solutions the series (12) can 
be approximated to p ~ q-\- /3\/{z 2 — z 2 ), the errors thus introduced being not more than those in- 
herent in all asymptotic expansions. This approximation, however, ignores the subtle difference 
between the polarization ratios of the upgoing and downcoming waves. 

Substituting this asymptotic form of p(z) in (11) and simplifying, it can be shown that the gen- 
eral type of magneto-ionic propagation (excluding a small region around = 77-/2) is characterized 
by the wave equation 



d 2 E x 

- + X 2 

dz 2 « 



Mg+j (4A/Kg) 
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ipEjp 
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where 



K 2 = ^7T 2 hya-qQ)l\ q = ±l 



M q 



_ iK q l 2 
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hi sin 2 6 
2(1 -qQ) 



1 

4' 



KH*h 2 sin 2 

A a = \ 3 



h 2 (l + sin 2 6) - qh 2 h 3 cos 6 - hi + 



qh\ cos 3 6 



2d -m l 

Equation (14) being similar to the transverse e-mode wave equation discussed in part I, the 
same procedure may be used for evaluating the phase factor Re^. Carrying out the simplifica- 
tions and neglecting terms of the order of h\ and smaller, it can be shown that 
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FIGURE 1. Reflection coefficient as a function of frequency. 
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FIGURE 2. Absorption coefficient as a function of frequency. 
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Figure 3. Time delay as a function of frequency. 
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where R 2 is the reflection coefficient of the g-mode for a friction-free layer, x q is its total absorption 
in the layer and r Q is the observed time delay in the CRO screen. The o-mode is q — + 1, while 
the e-mode is q = —\. Variations of/ 9 , x g /(gi>/2), and y Q /g with a>/a> m for 0=150°, 135°, and 120° 
have been plotted in figures 1, 2, and 3. The o-modes are shown by solid curves, and the e-modes 
by dotted curves. The two modes exchange their roles when 6 is replaced by (6 — tt/2). 

4. Polarization of the Downcoming Waves 

The electromagnetic (em) waves inside a magneto-ionic medium are not transverse with respect 
to the electric vector; hence, it has become customary to study the polarization of the waves with 
respect to the magnetic vector, which for vertical propagation is always confined to the horizontal 
plane. Since E • H = even for an inhomogeneous magneto-ionic medium, 

H x IH y =-EylE x = ip=^e i{ ^ ) , (16) 

where p is given by (12) and (13). Equation (16) shows that the vector whose components are H x 
and H y traces an ellipse in the %y-plane, the ratio of the major to the minor axes being given by ^, 
whereas the inclination ¥ (tilt angle) of the major axis to the magnetic meridian (xz-plane) is given by 

„_ 2$ sin \b 
t» 2¥ = 7r ^f (17) 

While evaluating the expression for the complex phase factor Re~ i<b , retention of terms up 
to /3i/(z 2 — Zq) in (12) is sufficient to predict results within the limits of experimental accuracy. In 
the polarization studies, however, terms up to /32l(z 2 — z 2 ) 2 need be retained. 

Integration of (10) when p = q+ 2 _^_ 2 + 2 __ 2 2 , 2 gives the dominant term as 

z z () (z z ) 

(18) 



where 



E 2_9< 


>*Bexp[g„z 2 /2] 


y/(dpjdz) 
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Since z is measured from the tip of the layer, the downcoming wave is given by that root of 
/^2 for which Im g ^ 0. Hence the limiting polarization p t for the downcoming wave is given by 



IR2 






(Z 2 -z§) '(l 2 -zl) 2 



J-T^+fva-w^o, 



(19) 



Separating the real and imaginary parts of (19) and neglecting terms containing v 2 /o) 2 or (tv/co, 
it can be shown that 



(20) 



* V ■ '- ■ 2 q 7 




-Zya41+ya>c(l+4y*)lq] 




v " m q 


y = co///co, yL = ycos0, 


yr = y sin 0, 



(l + 2y$)V(l-yf)(l + qy,.) , (21) 

(o c = sin 2 0/2 cos 9. 

The expression for c j#„ the ratio of the axes of the polarization ellipse, is the same as that 
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deduced by the ray optics, where electron density is assumed to be constant. Expression (21) 
for the tilt angle ^ contains two terms, W \ q a v, the collision frequency, and ^^ cr, i.e. \j{l\ m ). 
The second term, which depends on the characteristics / and a) m of the layer, is absent in the 
equivalent expression deduced from the ray optics. For ionospheric layers, the second term y i r 2q 
is greater than the first term ^i g , i.e., for these layers the effect of the nonhomogenity on the tilt 
angle ^ is more pronounced than the effect of the collision frequency v on it. If ^ and ^ e rep- 
resent the tilt angles for the o- and e-modes respectively, then 

Vo+^-^f-^r (22) 

o)fj I cos 6 

Experimental results on these tilt angles can thus be used to predict the value of /, the half-width 
of the layer. 

5. Conclusions 

The propagation of em waves in a magneto-ionic medium with a parabolic electron distribu- 
tion has been discussed in detail, and expressions for the observables have been derived. The 
coupled wave equations for the general station do not give rise to any third echo as expected by 
Rydbeck [1950]. 1 For other types of layers, or with collision damping as a function of the vertical 
height, the nature of the solutions will be different from those given in the text of the paper, though 
the basic predictions will not be affected to any marked extent. Expressions for the reflection 
coefficient and the ratio of the axes of the polarization ellipse agree, within limits of experimental 
errors, with those deduced from the ray optics. The equivalent height-frequency curve for a 
parabolic layer as deduced from the ray method would be an ellipse. It is not so from the wave 
optical standpoint. Similarly, the expressions for the tilt angle, as deduced by the two methods, 
differ radically from each other. The sharing of the energy between the two modes on entering 
the layer is, as discussed in the body of the paper, determined by the characteristics of the station. 
This method of evaluating the asymptotic solutions, and thereby the expressions for the observ- 
ables, can be utilized for any other type of electron layer as well. 

6. Appendix: List of Symbols and Abbreviations 

Field vectors E , H. 

Permittivity of free space eo- 

Permeability of free space /jl 

Steady magnetic field H e 

Collision frequency v 

Velocity of light c 

Charge and mass of electrons — e, m 

Electron density TV 

Maximum electron density N m 

Half- width of the layer / 

Propagation angle 6 

Magnetic dip angle 8 

Vertical height of the bottom of the layer hb 

Observed time delay on the CRO screen r q , To, r e 

Reflections coefficient for friction-free layer R 2 . R^ R 2 e 

Total absorption by the layer Xq-> Xo> Xe 

Ratio of the axes of the polarization ellipse &l q 

Tilt angle of the polarization ellipse ^ q 

1 O. E. H. Rydbeck (1950), Natl. Acad. Sci. (India) 4, 329. 
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Abbreviations: 
un= enoHe/m, 

&> 2 , = N m e 2 /me 

= 47T 2 C 2 /X m , 

h\ = l 2 l(l-iZ)Kl, 
_ (z*-a*)(z*-b*) 

(z 2 -^ 2 ) 
/^ 2 = 477- 2 M/(l-(? 2 )/ 4 , 



y=(t) H l(o, 
Z — v/to, 

hl = {l-iZWIa>%, 



yr = y sin 0, Jl = J cos 0, 
(o r = sin 2 0/2 cos 0, Q = h 3 cos 0/A 2 

M = ^/(l-^)w 2 m , 



H 



(z 2 -g 2 )(z 2 -6 2 ) 
(z 2 -z 2 ) ' 
/Y 2 sin 2 0'1 

1 _ /,2 -|_ _^ 72 

H 2 ^ l_Q-2 l ' 



L _^ (z 2 -a 2 )(z 2 -Z 2 ) 



U 2 -zf.) 



CO c 



A 2 cos 2 



a 2 = (l-Al)/ 2 , af = l-Af*- 2 -^ — -+A 2 A :i A/sin 2 + 



(? 2 cos 2 



/ 2 , 



^(l-Af + Ai)/ 2 , &§ = 



[l-Af + ^^-A 2 A^ 



cos 2 . , / . 9ni (? 2 cos 2 
"2 Ma y sin 2 + 



l\ 



M 1 , 2 =yK 2 + 6f,2-4), M = ^[a 2 + / 2 -z 2 ], 



A 1>2 =^(z8-a?, 2 )(z8-6f, 2 ), A = ^(zg-a?)(4-/ 2 ). 



£=iKz 2 , £o=iKzl 2/x=M, + M 2 , 



M Q 4 £ 



1 19,11 /I ^3 Sin S 

1 — hi + g/i 2 /i3 cos — — — 

2(1 — q(J)_ 



_ K 2 l 4 h 2 sin 2 



A§( 1 + sin 2 0) - gA 2 A 3 cos - A 2 + 



qh% cos 3 0] 



2(1 -Q 2 ) 
Kl=Wh\kl-qQ)l\ 

[1 — (co 2 jcx) 2 m + cof, sin 2 0/co 2 w )][l — qtOH cos 0/a>] 



/«=' 



y q =g 



V [1 — gco// cos 0/a>] 
l-|+|ln d)j/V[l- 9 a» w cos 0/a»], 



X q = {gvj2) 



1-|W|| 



[1 — gco« cos 0/w] 3 ' 2 , 



0, 2g/3 2 " 

2 9 ( 2 g-af)-y + — 
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q h h 72 • 2 a/9 a / ft>ft)/i\ / 2 sin 2 fl 

pi = h2h 3 l 2 sin 2 0/2 cos 0= — r-l -> 

\ (oli ) 2 cos 



{i-T^} + fVU-«i + ,« 



^i, =— Zyco, 



l + ^(l + 4y?) 



#», 



*2 9 = - co'yMd + 2y|)V(l-yf)(l + m )]/^ (J , 



%={i 



%=\ q+y^c+Yy 2 * * 
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